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THEORETICAL PERFORMANCE OF LIQUID HYDRAZINE AND LIQUID FLUORINE 

AS A ROCKET PROPELLANT 
By Sanford Gordon and Vearl N. Huff 



SUMMARY 

Theoretical values of performance parameters for liquid hydrazine 
and liquid fluorine as a rocket propellant were calculated on .the 
assumption of equilibrium composition during the expansion process for 
four chamber pressures and a wide range of fuel-oxidant and expansion 
ratios. The parameters included were specific impulse, combustion- 
chamber temperature, nozzle-exit temperature, equilibrium composition, 
mean molecular weight, characteristic velocity, coefficient of thrust, 
ratio of nozzle-exit area to throat area, specific heat at constant 
pressure, coefficient of viscosity, and coefficient of thermal conduc- 
tivity. A correlation is presented for the effect of chamber pressure 
on several of the parameters . 

The maximum value of specific impulse was 315.3 pound-seconds per 
pound for a chamber pressure of 300 pounds per square inch absolute 
(20.41 atm) and an exit pressure of 1 atmosphere. 


INTRODUCTION 

Liquid hydrazine and liquid fluorine are of interest as a rocket 
propellant because of high performance. -Extensive data exist in the 
literature on their availability, cost, and physical, chemical, and 
handling properties. 

The perfor ma nce of liquid hydrazine and liquid fluorine has been 
reported in the literature by a number of organizations such as Jet 
Propulsion Lab., C.I.T.; Douglas Aircraft Co., Inc., Project Rand; and 
North American Aviation, Inc. Additional performance calculations for 
this propellant covering a wide range of conditions were made at the 
NACA Lewis laboratory as part of a series of calculations on propellants 
containing the chemical elements hydrogen, fluorine, and nitrogen 
(refs. 1 and 2) to provide a comparison with the performance of other 
propellants based on the same thermodynamic data and computed to the 
same degree of accuracy, and to provide several parameters not pre- 
viously published. 




Data were calculated on the basis of equilibrium composition dur- 
ing expansion for four chamber pressures and a wide range of fuel- 
oxidant and expansion ratios. The performance parameters included are 
specific impulse, combustion-chamber temperature, nozzle-exit tempera- 
ture, equilibrium composition, mean molecular weight, characteristic 
velocity, coefficient of thrust, ratio of nozzle-exit area to throat 
area, specific heat at constant pressure, coefficient of viscosity, 
and coefficient of thermal conductivity. 


A correlation was made which permits the determination of specific 
impulse, characteristic velocity, and ratio of nozzle-exit area to 
throat area for a wide range of chamber pressures. Equations are given 
that permit the calculation of specific impulse for nonisentropic 
expansion or for change in heat content of the combustion gases using 
the originally calculated data. 

Several additional calculations were made assuming frozen compo- 
sition so that data based on the assumptions of equilibrium and frozen 
composition during the expansion process could be compared. 


SYMBOLS 

The following symbols are used. in this report: 

A number of equivalent formulas (function of pressure and molecular 
weight; see ref. 3) 

a local velocity of sound, ft/sec 

Cp coefficient of thrust 

C° molar specific heat at constant pressure, cal/(mole) (°K) 

Cp specific heat at constant pressure, cal/(g)(°K) 

c v specific heat at constant volume, cal/(g)(°K) 

c* characteristic velocity, ft/sec 




g 


=? 

h 

I 

k 

M 

n 

P 

P 

<1 

R 

r 

S 

s 

T 

w 

Y A 

Y i 

T s 


acceleration due to gravity, 32.174 ft/sec 2 

sum of sensible enthalpy and chemical energy, cal/mole 

sum of sensible enthalpy and chemical energy per unit weight, 

-hs — cal/g 

specific impulse, lb-sec/lb 

coefficient of thermal conductivity, cal/(sec) (cm) (°K) 

molecular weight, g/mole 

number of moles j exponent 

pressure 

partial pressure 

heat, cal/g 

universal gas constant (consistent units) 

equivalence ratio, ratio of number of fluorine atoms to hydrogen 
atoms 

nozzle area, sq ft 
entropy, cal/(g)(°K) 
temperature, °K 
rate of flow, lb/sec 




H coefficient of viscosity, g/(cm)(sec) ■ poise 

p density, g/cu cm 

Subscripts: 

c combustion chamber 

e nozzle exit 

frozen composition assumed frozen 
i product of combustion 

max maximum 

o initial value 

P constant pressure 

s constant entropy 

t nozzle throat 

x any point in nozzle 


CALCULATION OF PERFORMANCE DATA 

Calculations of the performance data were made with a Bell computer 
and an IBM Card -Programmed Electronic Calculator as described in refer- 
ence 1. The assumptions, thermodynamic data, and transport properties 
used for the calculations are the same as those of reference 1. 

The products of combustion were assumed to be ideal gases and 
included the following substances: hydrogen fluoride HF, hydrogen 

H 2 , nitrogen Ng, fluorine F 2 , atomic fluorine F, atomic hydrogen H, 
and atomic nitrogen N. The dissociation energy of F 2 was taken to 
be 35.6 kilocalories per mole (ref. 4). Physical and thermochemical 
properties of the propellants were taken from references 3 to 6 and are 
given in table I. 

Procedure for combustion conditions . - The following parameters 
were computed for ten equivalence ratios for a chamber pressure of 
300 pounds per square inch absolute and for five equivalence ratios for 
chamber pressures of 150, 600, and 1200 pounds per square inch absolute: 
combustion temperature, equilibrium composition, enthalpy, mean 
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molecular weight, derivative of the logarithm of pressure with respect 
to the logarithm of density at constant entropy r s > specific heat at 
constant pressure, coefficient of viscosity, coefficient of thermal 
conductivity, and entropy of the combustion products. 

Procedure for exit conditions . - Equilibrium composition, mean 
molecular weight, pressure, derivative of the logarithm of pressure with 
respect to the logarithm of density at constant entropy r s , enthalpy 
of the products of combustion, specific heat at constant pressure, 
coefficient of viscosity, and coefficient of thermal conductivity were 
computed for each equivalence ratio by assuming isentropic expansion 
for four assigned exit temperatures selected to cover the exit pressure 
range from the nozzle-throat pressure to about 0.1 to 0.3 atmosphere. 

Interpolation. - Parameters for pressures at and near the nozzle 
throat, for pressures of 2, 1, 0.5, 0.25, 0.125 atmosphere, and for 
pressures corresponding to altitudes of 10,000, 20,000, 30,000, 40,000, 
and 50,000 feet were interpolated by means of cubic equations between 
each pair of the assigned exit temperatures. The functions and their 
first derivatives used in the interpolations are described in refer- 
ence 1. 


The errors due to interpolation were checked for several cases. 
The values presented for all performance parameters appear to be cor- 
rectly interpolated or in error at most by two or three units in the 
last place tabulated. 


Formulas. - The formulas used in computing the various parameters 
are given in reference 1 and are summarized as follows: 


Specific impulse, lb-sec/lb 

I = 294.98' 


I 1000 


( 1 ) 


Throat area per unit flow rate, (sq ft)(sec)/lb (pressure in atm) 


S t /v = 


1.3144 T t 
a 


Characteristic velocity, ft/sec 


c* = g P c S t /w = 32.174 P c S t /w 


( 2 ) 


( 3 ) 







Uarc-Sd 



Coefficient of thrust 


C F = Ig/c* = 32.174 i/c* ( 4 ) 

Nozzle-exit area per unit flow rate, (sq ft)(sec)/lb (pressure in atm) 

„ ; 0.040853 T P 

S e /w = 


P e M e I 


( 5 ) 


( 6 ) 


Ratio of nozzle-exit area to throat area 

s /s = Se/w 
e t S^./w 

Specific heat at constant pressure, cal/(g)(°K) 

C P = [ T ^ n i (C P 5 i + 2>i(H?)i Yi -ZoiC^i Y A ] (7) 

Derivative of the logarithm of pressure with respect to the logarithm 
of density at constant entropy 

I>i Di 


^s " 


Pfo A -l) 


( 8 ) 


Coefficient of viscosity, poise 


PM 


( 9 ) 


P V' p i 

^(,i/Mi) 

Coefficient of thermal conductivity, cal/(sec) (cm) (°K) 

* 3 4* + f !) (10) 

When composition is assumed to be frozen, the partial derivatives 
and Y A in equation (7) are equal to zero and the partial deriva- 
tives Di and D A in equation (8) are equal to — p ^° Z9 P . There- 
fore, equations (7) and (8) become: 

Specific Q heat at constant pressure assuming frozen composition, 

Cal/(g)(K) ■ SMeJh 

( c p^frozen 37 (-*■-*-) 


nM 
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Derivative of the logarithm of pressure with respect to the logarithm 
of density at constant entropy assuming frozen composition 


(O 


( c p) frozen 


s' frozen 




v p'frozen 


r/m 



frozen 


( 12 ) 


The values of viscosity and thermal conductivity for mixtures of 
combustion gases calculated by means of equations (9) and (10) are 
only approximate. When more reliable transport properties for the var- 
ious products of combustion become available., a more rigorous procedure 
for computing the properties of mixtures may also be justified. 


THEORETICAL PERFORMANCE DATA 

The calculated values of the performance parameters specific 
impulse , temperature, mean molecular weight, characteristic velocity, 
coefficient of thrust, and ratio of nozzle-exit area to throat area are 
given in tables II to IV. For convenience, the parameters are tabu- 
lated as functions of altitude in table II and as functions ‘of expansion 
ratio in table III. As an aid to engine design, the values of the par- 
ameters within the rocket nozzle for 80, 90, 100, 110, and 120 percent 
of the throat pressure are tabulated in table IV. Equilibrium compo- 
sition, y q > specific heat at constant pressure, coefficient of viscos- 
ity, coefficient of thermal conductivity, and mean molecular weight in 
the combustion chamber and at assigned exit temperatures are given in 
table V. The mole fraction of Fg was not tabulated in table V, since 
it was always less than 0.00003 except at 1600° K for the equivalence 
ratio r » 1.2, where it was 0.00007, 0.00015, and 0.00030 for the 
chamber pressures of 300, 600, and 1200 pounds per square inch absolute, 
respectively. 

Parameters . - The parameters are plotted in figures 1 to 13. 

Curves of specific impulse for various altitudes are shown in figure 1 
plotted against weight percent fuel. The maximum value of specific 
impulse for the sea-level curve in figure l(b) is 315.3 pound-seconds 
per pound at 31.9 percent fuel by weight. 

The maximum values of specific impulse and the weight percentages 
at which they occur were obtained by numerical differentiation of the 
calculated values and are shown in figure 2 as a function of altitude. 
The maximum specific impulse in figure 2(b) increases 22 percent for 
a change in altitude from sea level to 50,000 feet. 



ui ^ t 



Combustion-chamber temperature and nozzle-exi't temperature for 
various altitudes are presented in figure 3 as functions of weight per- 
cent fuel. The maximum combustion temperature in figure 3(b) is 
4411° K at 27.7 percent fuel by weight. The maximums of the exit tem- 
perature curves occur near the stoichiometric ratio. 

Characteristic velocity and coefficient of thrust are plotted in 
figure 4 and the ratio of the area at the nozzle exit to the area at 
the throat is plotted in figure 5 against weight percent fuel. 


Curves of mean molecular weight in the combustion chamber and 
nozzle exit are plotted against weight percent fuel in figure 6. 

Curves of specific heat at constant pressure, coefficient of vis- 
cosity, and coefficient of thermal conductivity for six pressures are 
plotted in figures 7 to 9 as functions of weight percent fuel. 


Chamber pressure effect . - The values of the parameters I, c*, 
and Sg/S-t- (or the logarithms of these parameters) are very nearly lin- 
ear with the logarithm of chamber pressure for a fixed equivalence 
ratio and expansion ratio (for example, see fig. 10) . It was therefore 
possible to correlate these data for chamber pressure effects according 
to the following equations to the accuracies indicated: 


1 = I 300 (§o) ( ±0 - 06 Percent) 

c * = c 300 (350) ( ±0 * 07 Vercent) 

S e /S t = (S e /S t ) 3 00 (§ 5 ) ( ±0 - 3 Percent) 


(13) 

(14) 

(15) 


where 


T e* 
■*■300 ^ c 300 


, and C^ e /®t^300 are va -*- ues °f these parameters 


at 300 pounds per square inch absolute, P c is in pounds per square inch 
absolute, and n is an exponent given as a function of fuel-oxidant 
and expansion ratios for each parameter. The values of n shown in 
figures 11 to 13 were computed from the values of the parameters for 
the various chamber pressures. However, values of n for specific 
impulse can also be calculated from data for a single chamber pressure 
by means of equation (58) derived in the appendix. All the data given 
in this report for I, c*, and S e /S+ are correlated in figures 11 
to 13. 


As an illustration of the use of these curves, suppose it is 
desired to obtain the values of I, c*, and Sg/S-t for a combustion 



pressure of 800 pounds per square inch absolute and an expansion ratio 
of 150 at 30 percent fuel by weight. The values read from figures 11 
to 13 are: 


l 300 


= 380 


c 300 “ 7090 

»A - 17 - 3 


n = 0.0072 
n = 0.0141 
n = -0.042 


Using these values in equations (13) to (15) gives Iqoo = 383 > 
c 800 = 7 ^ 89 > an< ^ ( 8 e/ 8 t^800 = 16.6 


Change in heat content of propellant gases . - The results pre- 
sented in this report are computed for adiabatic combustion with pro- 
pellants at the initial temperature indicated. For a change in heat 
content of the combustion gases in the combustion chamber (which would 
result, for example, from heat loss in the combustion-chamber or intro- 
duction of the propellants at a temperature other than the one indi- 
cated) , the corresponding specific impulse assuming isentropic expansion 
and the same combustion and exit pressures as in the initial calcula- 
tions may be obtained from the following equation: 

I 2 = I 2 + B Ah c + C (Ah c ) 2 (16) 


where Ah c is the change in heat content of the combustion gases in 
the combustion chamber, 


B = 87.0132 



C = 




and the subscript o indicates the values of the parameters before the 
change is made. The derivation of this equation is given in the 
appendix. 

For example, assume that the enthalpy of the propellant gases fol- 
lowing combustion is 100 calories per gram less than the initial value 
of h c used in these calculations for the stoichiometric mixture ratio 
at a chamber pressure of 300 pounds per square inch absolute 
(Ah c = - 100 ) . For an expansion ratio of 20.41 (P e = 1 atm), the values 
of I and T obtained from either table III (b) or figures l(b) and 
3(b) and the values of Cp obtained from figure 7(b) are as follows: 
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I 0 = 314.3 lb-sec/lb 
T= 4394° K 

) U 

T = 3239° K 
e,o 

(c p) C ,o = 1,958 cal /(g)(°K) 
( c p) e ,o = 1>122 ca l/(g)(°K) 
These values yield the following: 

I 2 = 98,784 


B = 22.87 
C = -0.00278 


By equation 16, 

I 2 = 98,784 + 22.87 (-100) + (-0.00278) (10,000) 


= 98,784 - 2,287 - 28 =,96,469 
I = 310.6 lb-sec/lb 

Equation (16) was applied to a specific case having a Ah c of 
250 calories per gram and resulted in an error of less than 0.1 pound- 
second per pound in specific impulse. For other cases, the relative 
magnitudes of the successive terms in the series may be used to esti- 
mate the approximate magnitude of the error. 

Nonisentropic expansion . - For design of rocket engines and for 
correction of experimental results for heat loss from the nozzle of 
experimental engines, a need exists for performance data assuming non- 
isentropic expansion. In general, the performance assuming nonisen- 
tropic expansion through the nozzle cannot be computed until sufficient 
information is available to determine the heat rejection from the nozzle 
q and the entropy change during expansion As. When suitable informa- 
tion is available the computation of performance is possible with the 
aid of an enthalpy-entropy diagram or by use of an equation such as the 
following one using data calculated for isentropic expansion: 


I 2 = I 2 - 87.0132 + T e q (As) + 


This equation is derived in the appendix. 


(c ) 

L p'ej 


(As)‘ 


(17) 



The use of equation (17) is illustrated by imposing an arbitrary 
relation between q and As by assuming that all the heat q is 
removed from the nozzle throat at constant pressure. For the stoichio- 
metric mixture ratio at a chamber pressure of 300 pounds per square 
inch absolute and an expansion ratio of 20.41, the values of T^. 
obtained from table IV(b) and of ( c p)t obtained from figure 7(b) are 
as follows : 

T = 4163° K 

"C 

(c p ) t = 1.860 cal/(g)(°K) 

If q = 100 calories per gram, then by equation (44) of the appendix. 

As = -0.02402 - 0.00016 = -0.02418 cal/(g)(°K) 

This value of As, the value assumed for q, and the values of I Q , 

Te,o> an 8 ( c p) e o given in the previous example are used in equa- 
tion (17) to give the following result: 

I 2 = 98,784 - 87.0132 jlOO + 3239 (-0.02418) + 2887 

= 98,784 - 87.0132 (100 - 78.32 + 0.84) = 96,824 

I = 311.2 lb-sec/lb 

Frozen composition . - In order to compare data based on the assump- 
tions of equilibrium and frozen composition during the expansion pro- 
cess, several additional calculations were made assuming frozen com- 
position. These are presented in table VI together with corresponding 
equilibrium data for the stoichiometric equivalence ratio and for two 
expansion ratios. The percentage differences in these parameters for 
frozen and equilibrium composition are considerably higher for the 
expansion ratio of 163.3 than for the expansion ratio of 20.41. 

For a combustion pressure of 300 pounds per square inch absolute 
and an exit pressure of 1 atmosphere, the values of maximum specific 
impulse are 315.3 pound-seconds per pound at 31.9 percent fuel by weight 
for equilibrium composition during expansion and 292.2 pound-seconds 
per pound at 32.5 percent fuel by weight for frozen composition during 
expansion . 




SUMMARY OF RESULTS 


Theoretical calculations of the performance parameters of liquid 
hydrazine with liquid fluorine were made for a wide range of chamber 
pressures, expansion ratios, and fuel-oxidant ratios and yielded the 
following results: 

1. For a combustion -chamber pressure of 300 pounds per square inch 
absolute (20.41 atm) and an exit pressure of 1 atmosphere, the maximum 
specific impulse was 315.3 pound-seconds per pound at 31.9 percent fuel 
by weight for equilibrium composition during expansion and 292.2 pound- 
seconds per pound at 32.5 percent fuel by weight for frozen composition 
during expansion. 

2. The maximum combustion temperatures calculated for chamber pres- 
sures of 150, 300, 600, and 1200 pounds per square inch absolute were 
4260°, 4411°, 4539°, and 4687° K, respectively. 

3. The percentage differences in the parameters specific impulse, 
coefficient of thrust, ratio of nozzle-exit area to throat area, exit 
temperature, and exit mean molecular weight due to the assumption of 
frozen or equilibrium composition during expansion are considerably 
higher for the expansion ratio of 163.3 than for the expansion ratio 
of 20.41. 

4. All the data presented in this report for specific impulse, 
characteristic velocity, and ratio of nozzle-exit area to throat area 
have been correlated for chamber-pressure effect in figures 11, 12, and 
13, respectively. 

5. An equation is given for computing the change in specific 
impulse resulting from a change in the heat content of the combustion 
gases . 


6. An equation is given for computing the value of specific impulse 
for a nonisentropic expansion from data for isentropic expansion. 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, April 9, 1953 


APPENDIX - DERIVATIONS OF EQUATIONS FOR NONIDEAL PERFORMANCE 
AND FOR EXPONENT IN PRESSURE CORRELATION 

All the data presented in this report were computed on the assump- 
tions of adiabatic combustion and isentropic expansion with propellants 
at the initial temperature indicated. Figure 14 illustrates these 
assumptions on an enthalpy- entropy diagram. In practice, however, 
these ideal assumptions are not realized. Equations are derived, there- 
fore, which permit using the data in this report to compute specific 
impulse for a change in the enthalpy of the propellant gases (eq. (29), 
see fig. 15), or for nonisentropic expansion (eq. (37), see fig. 16); 
or to compute the change in entropy for a change in enthalpy along a 
constant pressure line (eq. (44), see fig. 17). 

Correction of specific impulse for a change in enthalpy of propel- 
lant gases. - As may be seen from figure 15, h e is a function of he 
for isentropic expansion from a fixed chamber pressure to a fixed exit 
pressure , 



h e > fog 



(19) 

Expanding this function 

in a Taylor's 

series 

and neglecting all 

derivatives higher than the 

second yield 







Ah 2 

(20) 

- f C h c,o> + 

f ' (h ) Ah + 
' c,o' c 

f"(n C) 

\ \ c / 
o> — 2 — 

For constant pressure , 






dh = T ds 



(21) 


Oj\Oj 

tHi-3 

!! 



(22) 


For isentropic expansion , 

ds„ , % 

s a s or - — = 1 (23) 

e c ds c 


The relations in equations (21) , (22), and (23) are used to obtain 


f’M 


T e ds e _ 

dh c r ^c ( ^ s c T c 


(24) 
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f"(hc) = 


1 l^A m 
T c V dh J + 6 


T c dh e dh c + T e 


m 2 
X C 


Al \ ® c m l m 1 (V\ _ ^e_ _1 

\T C 2 / dh c T c ( c p)e \ T c / T c 2 ( c p 

r_i lj 

_^ c p^e ^ c p^cJ 


p^c 


(25) 


Substituting equations (19), (24), and (25) into equation (20) gives 

... ^2 


h e = 


^ + 

1 1 

/O \T c 2 / 0 

( C p^e ( c p ) c 


(26) 


From equation (l), 

I 2 = 87.0132 (h c - he) (27) 

By definition, 

= h c ' h c,o ( 28 ) 

Substituting equations (26) and (28) into equation (27) results in 


I 2 = 87 


•0132 l^he + h c , 0 - ^e,o ” ^T c ) 0 


\ m 2/ 

v X C / 0 


^ c p^e xy; 


(Ah c )‘ 


By use of equation (27) this may be expressed 
2 t 2 . e ., U T e\ a t_ . 87.0132 / T e\ r 


1=1+ 87.0132 
o 


(*-a 


Ah c + 2 


C Jo 


T^Tc'T^l 


CtocV 


Equation (29) is identical to equation (16) in the text. 


BaBiisgssssa ^ 


mstasstsmr‘ 


(29) 
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Correction of specific impulse for nonisentropic expansion . - As 
may be seen from figure 16 , h e is a function of s g for nonisentropic 
expansion, since h and s. must lie on the line P . 

“ G 0 

h e = f ( s e ) ( 3 °) 

Expanding this function in a Taylor's series and neglecting all deriv- 
atives higher than the second give 

h = f(s o J + f ' (s Q J As + f"(s ) ( 31 ) 

e x e,o' ' e,o' ' e,o' 2 ' ' 

For constant pressure, 

f ,( s ) = || = T (32) 

f,, ( s > - = S - f . < 33 ) 

ds 2 as C P 

Substituting equations (32) and (33) into equation (3l) yields 


h e " h e,o + ( T e,o) As + 


m —1 

T e 

( c p)e o 


(As) 2 


When equation (27) is modified to account for heat q lost during 
expansion, the result is 

I 2 = 87.0132 (h c - h e - q) (; 

Taking h c ^ Q = h c and substituting equation (34) into equation (35) 


I 6 = 87.0132 <h - h - (T ) As 
| c,o e,o v e,o' 


T e 

_( c ph_o 


^--4 (36) 


Equation (27) is used to express this as 


I 2 = I 2 - 87.0132 <q + (T ) As 
o 1 e,o 


+ T T e 1 

Try 

' p'e c 


(As) 2 


Equation (37) is identical to equation (17) in the text. 

Change of entropy for change of heat along a constant pressure 
line. - As may be seen from figure 17, for a change in enthalpy along 
any line of constant pressure P x , 
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s x “ f < h x) 


Expanding this function in a Taylor's series and neglecting all deriva- 
tions higher than the second result in 

s x - + f, ( h x,o) ^ + f'OSc.o) < 39 ’ 

For constant pressure , 


f"(h) 


- & 


i 

T 


O d f — | 

dfs _ It/ -1 dT 

dh »p2 dh 



(40) 

(41) 


Equations (38) , (40) , and (4l) are substituted into equation (39) to 
give 



(42) 


When the change in enthalpy at constant pressure P x is due to 
removal of heat q from the system, 

q = -Ah (43) 

(Addition of heat would correspond to negative q.) 

With the aid of equation (43) , equation (42) becomes 



An example is given in the text of the use of this equation when 
a nonisentropic expansion is approximated by assuming that the heat is 
removed at an average pressure. 

Equation (44) may also be used to obtain additional points on the 
constant-pressure lines of an enthalpy-entropy diagram. 

Calculation of chamber pressure correlation exponents . - In the 
general notation, equation (13) may be written 
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I * I 


n 


0 VP, 


c .o 


(45) 


Values of n to be used in equation (45) for obtaining specific 
impulse at any chamber pressure P c from the value of specific impulse 
at an initial chamber pressure P c 0 for the same expansion ratio and 
equivalence ratio may be obtained as follows: From equations (27) and 

( 45 ), 

\2n 


I 2 = 87.0132 (h c -h e ) = I 2 


Differentiating equation (46) with respect to P,, gives 

(p„ 


87 - 0132 (a£ - S?) - 2n 


L c 
2n-l' 


2n 


c ,o 


Since h c is constant, 


dh = 0 
c 


Substituting equations (46) and (48) into equation (47) yields 

(& h e \ 

- 87.0132 —5 = 2n i- 

VdP c ; P c 


(46) 


(47) 


(48) 


or 



(49) 


Since entropy s may be expressed as a function of P and h. 



Substituting equation (48) into equation (50) results in 
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ds < 3 te! c ap = (52) 

For isentropic expansion ds e = ds c ; therefore from equations (5l) and 
(52), 



For constant expansion ratio, 



Substituting equation (54) into equation (53) gives 



(53) 


(54) 


or 



Substituting into equation (55) the thermodynamic identity of 
equation (40) and the following identity 


give 



dhe _ fi_ l_\ 

dp c ~ " p c \ lA c ' M e ) 

Equation (57) is substituted into equation (49) to yield 


(56) 


( 57 ) 
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The following table presents values of n computed by means of 
equation (58) for the stoichiometric mixture ratio for a chamber pres- 
sure P c o of 300 pounds per square inch absolute. Also presented is 
a comparison of the values of specific impulse at other chamber pres- 
sures computed directly and computed by means of equation (45) from 
these values of n and 1300* 


p c/ p e 

n 

(from 
eq. (58)) 

P c, 

lb/sq in. abs 

I 

(direct com- 
putation) , 
lb-sec/lb 

I 

(from 

eq- (45)), 
lb-sec/lb 

Error 
in I, 

lb-sec/lb 

20.41 

0.01240 

150 

311.53 

311.58 

0.05 



300 

314.27 


— 



600 

316.92 

316.98 

.06 



1200 

319.45 

319.72 

.27 

81.65 

0.00966 

150 

361.27 

361.37 

0.10 



300 

363.80 


— 



600 

366.16 

366.24 

.08 



1200 

368.33 

368.70 

.37 

163.3 

0.00804 

150 

380.48 

380.59 

0.11 



300 

382.72 





600 

384.76 

384.86 




1200 

386.62 

387.01 



Equations similar to equation (58) may be derived for evaluating 
the exponent n for c* and S e /S t ; however, these equations are 
omitted inasmuch as they involve partial derivatives which have not 
been tabulated in this report. 
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TABLE II. - CALCULATED PERFORMANCE OF LIQUID HYDRAZINE WITH LIQUID FLUORINE AS FUNCTIONS OF ALTITUDE 
(a) Combustion-chamber pressure, 150 lb/sq in. absolute 
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^ased on F 2 density of 1.54 at -196° C and NgH^. density of 1.011 at 15° C. 

^Nozzle designed for exit pressure equal to ambient pressure corresponding to altitude Indicated. 



TABLE II. - CALCULATED PERFORMANCE OF LIQUID HYDRAZINE WITH LIQUID FLUORINE AS FUNCTIONS OF ALTITUDE - Continued 

(b) Combustion-chamber pressure, 300 lb/sq in. absolute 


NACA EM E53E12 



Specific 

impulse, 

I, 

lb-sec/1 

CD co co r- cd co 

co cd cd ft in in 
o ft to ft in co 
to ro to to to to 

CD cd O ft cd r- 
ft in o to in co 

ftCDftincor- 
to to to to to to 

tO ft CD CO ft CO 
ft CD ft CO CD ft 

ft cd ft in t~- co 
to to to to to to 

tO CD > o co 

in o ft co cd to 
ft to ft in r~- co 
to to to to to to 

ft O O CO ft CO 
ft CD to CO CD O 

ft cd ft in co co 
to to to to to to 

i - 

ft P 

r- o CD in CO CD 

ft cu co r- in ft 

r- co to co ft r-- 

ft CD r- O ft ft 

CD LO CD O CO CD 

ft p m 

ft co to CD ft CD 

CD CD L0 ft r- CD 

CD CD CO CD CD ft 

CD CD in CD co to 

ft co ft ft co ft 

ft C 73 Ph 

0) 0) u o 

o ft ft ^ 
o O O P 

ft ft in in to co 

ft ft in co co r- 

ft ft in co co t 

ft ft in to co r- 

ft ft LO CO CO r- 

r — l i — l rH rH rH rH 

ft ft ft ft ft ft 

ft ft ft rH ft ft 

ft ft ft ft ft ft 

ft ft ft ft ft ft 

of 

0 

ft ft ft 00 

CO CD tO CO 

CD CD ft in 

cd o in ft 

ft to co CD 


ft ft ft 
Cd -CO 


rl NP Bi O Cd\ 

P N ft Q) U CD <L> 

cd O * U S3 U CO 

CC £ <d cd p cd 


cor-cnintor» 
I CO CO O ft CD LO 


r-- to in co r- o 

CO O CO O) lO CM 


ft CD rH CD rH O 

cd h o r- co co 


go co o ft in to 

CO O 00 to I— I lO 


CD ft CD CD ft CO 

r- cd in o to r~ 


tO ft CO CO ft LO 


to in CO CO CD e'- 
en rH 


to LO CD to CD 1 to UO to CD to QO 


to ft co cd CD r- 
| HH 


I P 

o S3 

P <D U hO 0) 

Cd ft Cd ft S 

<D O I-H <1> 

S 6 73 3c 


in in rH CD CD CD 
CD to ft ft ft ft 


I LO CD co rH CD CO 
CD ft CD ft lO lO 


h in r- co to o 
in to r- co cd 


CDHHCDLOt^- 
r- CD O O ft ft I 


f-l - 
CD CD * 
(XU 0) w 
B 2 Eh O 
CD P 
EH Cd 


00 to CD CD CO CD 
rH r-- CD r- CD o 

CD CO ft ft CD f- 
CD CD CD CD rH H 


in in ft ft to ft 
r- 00 in CD CD r- 
ft cn t-- ft CD cd 
to CD CD CD CD rH 


cd cd m in co to 
to CD ft CO in to 
CD o CD LO to 
to to CD CD CD CD 


to to H CD CD CD 
tO t"- CD CD CO ft 

h cd r- in to h 

to CD OJ CD CD CD 


O^Hoor^r- 
1 CO CD H o co CO 
CD D- CO ft ft CD 
CD CD CD CD CD rH 


O >i 

I HP O 
O P ft - <D 

cd to o* ra 

JhH o o\ 
td U rH P 

S3 <D CD ft 

O P > 


1 p 

C 8 U bb o 

cd ft cd ft S 


CO ft CO CD CD 

r- cd co in ft 


CO ft CO CD CD 
r— cd co in ft 


co ft oo CD CD 
r- cd to in ft 

CO in CD CO rH 


CO ft 00 CD CD 

r>- cd to in ft 

CO in CD co rH 


CO ft CO CD CD 
r- cd co in ft 
CO in CD co rH 


to Ph cd 
Q) 

u 

pH 

o CO ft CD ft ft 
rH 

O CO ft CD i — 1 ft 
rH 

O co ft CD ft ft 

rH 

o co ft (Sj ft h 

ft 

O CO ft CD ft ft 
ft 

<D 

o o o o o o 

o o o o o o 

o o o o o o 

o o o o o o 

o o o o o o 

2 

o o o o o 

o o o o o 

o o o o o 

o o o o o 

o o o o o 

p p 

o o o o o 

o o o o o 

o o o o o 

o o o o o 

o o o o o 

ft ft 






P 

o o o o o 

o o o o o 

o o o o o 

o o o o o 

o o o o o 

IV 

ft cd to ft in 

ft cd to ft in 

ft cd to ft m 

ft cd to ft m 

ft cd to ft in 



cd 

u 

CD - - 
ft 0) o W 
E U EHO 
CD 2 
EH P 


to 

ft 

ft 

CD 

CD 

CD 

ft 

CD 

in 

r- 

to 

ft 

to 

to 

CD 

ft 

ft 

ft 

■ ft 

ft 


cd 

1 p 


ft 

P c 

o 

ft 

S3 0) 

o 

CD 

bO O ft 


p. 

ft U <D 

co 

o 

cd cd 3 

CD 

u 

32 (XU 


Pn 




cd - 



> <d o 

CD 


ft O ft U 

• 


2 C ft 

ft 


o' cu cd 



W ft u 



to 

o 

ft 

to 

CD 

o 

to 

to 

to 

ft 

ft 

1 — 1 


RESTR-KSE 


^ased on F 2 density of 1.54 at -196° C and N 2 H 4 density of 1.011 at 15° C. 

b Nozzle designed for exit pressure equal to ambient pressure corresponding to altitude indicated. 
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^Based on F 2 density of 1.54 at -196° C and N 2 H 4: density of 1.011 at 15° C. 

’Nozzle designed for exit pressure equal to ambient pressure corresponding to altitude indicated. 






NACA EM E53E12 



'O 

<D 

§ 

-P 

c 

O 

o 


8 


M 


O 

CO 

s 

o 

H 

Eh 

O 

P 


PC 

0 

1 


§ 


p 

Eh 

H 

W 

s 

H 

N 

< 


s 


O 


w 

o 


s 

PC 

o 

Ph 

« 

w 

pL| 

p 

w 

Eh 

< 

P 

P 

O 


w 

p 

g 

Eh 




Based on F 2 density of 1.54 at -196° C and NgH^ density of 1.011 at 15° C. 

'Nozzle designed for exit pressure equal to ambient pre.esure corresponding to altitude indicated. 
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a Based on F 2 density of 1.54 at -196° C and density of 1.011 at 15° C. 

b Nozzle designed for exit pressure equal to ambient pressure corresponding to altitude Indicated. 



TABLE III. - CALCULATED PERFORMANCE OF LIQUID HYDRAZINE WITH LIQUID FLUORINE AS FUNCTIONS OF EXPANSION RATIO 
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Based on F 2 density of 1.54 at -196° C and N 2 H 4 density of 1.011 at 15' 
Nozzle designed for exit pressure equal to ambient pressure. 










MCA EM E53E12 




o 


PL, 

X! 

W 

o 

00 

s 

o 

H 

Eh 

O 

S 3 

!=> 

Ph 

CO 

< 

W 

13 

H 

« 

o 

3 


B 

H 

W TJ 
Eh 0) 
H d 

^ c 

*H 
W -P 

s c 

H O 
N O 

S 1 

>H O 
« H 

Q § 
H 


<u 

-p 

3 


'S 


O 1 

w 

\ 


O' 

o 

fO 


0 ) 

5 h 

3 

03 

OQ 

<U 

?H 

a 

Jh 

< 1 ) 

1 

.c 

0 

1 

c 

o 

T-{ 

03 

d 

rQ 

S 

o 

o 





TABLE III. - CALCULATED PERFORMANCE OF LIQUID HYDRAZINE WITH LIQUID FLUORINE AS FUNCTIONS OF EXPANSION 

RATIO - Continued 

(b) Combust ion- chamber pressure, 300 lb/sq in. absolute - Concluded 
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^Based on P 2 density of 1.54 at -196° C and N 2 H 4 density of 1.011 at 15 
^Nozzle designed for exit pressure equal to ambient pressure. 



TABLE III. - CALCULATED PERFORMANCE OF LIQUID HYDRAZINE WITH LIQUID FLUORINE AS FUNCTIONS OF EXPANSION 

RATIO - Concluded 
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m Fg density of 1.54 at -196° C and N 2 H 4 density of 1.011 at 15° 
designed for exit pressure equal to ambient pressure. 
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TABLE V. - PROPERTIES AND COMPOSITION IN COMBUSTION CHAMBER AND FOLLOWING AN ISENTROPIC EXPANSION TO ASSIGNED EXIT TEMPERATURES FOR 

LIQUID HYDRAZINE WITH LIQUID FLUORINE - Concluded 
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TABLE VI. - COMPARISON OF CALCULATED PERFORMANCE OF LIQUID HYDRAZINE WITH LIQUID 
FLUORINE ASSUMING EQUILIBRIUM AND FROZEN COMPOSITION DURING EXPANSION 
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Figure 1. - Continued. Theoretical specific impulse of liquid hydrazine with liquid fluorine. 
Isentropic expansion to altitude indicated assuming equilibrium composition. 
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Figure 1. - Continued- Theoretical specific impulse of liquid hydrazine with liquid 
fluorine. Isentropic expansion to altitude indicated assuming equilibrium 
composition. 










Maximum specific impulse, I max , lb -sec/lb 



propellant, percent by weight 






(b) Combustion-chamber pressure, 300 pounds per 
square inch absolute. 

Figure 2. - Continued. Maximum theoretical specific impulse 
and corresponding weight percent fuel in propellant of liquid 
hydrazine with liquid fluorine. Isentropic expansion to alti- 
tude indicated assuming equilibrium composition. 
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(d) Combustion -chamber pressure, 1200 pounds 
per square inch absolute. 1 

Figure 2. - Concluded. Maximum theoretical 
specific impulse and corresponding veight 
percent fuel in propellant of liquid 
hydrazine ylth liquid fluorine. Isen- 
tropic expansion to altitude indicated 
assuming equilibrium composition. 
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Figure 3. - Theoretical combust ion -chamber temperature and nozzle-exit temperature of 
liquid hydrazine with liquid fluorine. Isentropic expansion to altitude indicated 
assuming equilibrium composition. 
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(c) Combustion -chamber pressure, 600 pounds per square inch absolute. 

Figure 3. - Continued. Theoretical combustion-chamber temperature and nozzle -exit 
temperature of liquid hydrazine with liquid fluorine. Isentropic expansion to 
altitude indicated assuming equilibrium composition- 
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Figure 5. - Theoretical ratloi 
with liquid fluorine- Iseni 
rium composition- 
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(b) Combustion-chamber pressure, 300 pounds per square inch absolute. 

Figure 5. - Continued. Theoretical ratios of nozzle-exit area to throat area for liquid 
hydrazine with liquid fluorine. Isentropic expansion to altitude indicated assuming 
equilibrium composition. 
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Fuel in propellant, percent by weight 

(c) Combustion -chamber pressure;, 600 pounds per square inch absolute- 


Figure 5. - Continued. Theoretical ratios of nozzle -exit area to throat area for 
liquid hydrazine with liquid fluorine. Isentropic expansion to altitude indicated 
assuming equilibrium composition- 
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(a) Combust ion -chamber pressure, 150 pounds per square inch absolute. 

Figure 6. - Theoretical mean molecular weight in combustion chamber and at nozzle 
exit for liquid hydrazine with liquid fluorine. Ipentropic expansion to altitude 
indicated assuming equilibrium composition. 











Fuel in propellant, percent by weight 

(c) Combustion-chamber pressure, 600 pounds per square inch absolute. 

Figure 6. - Continued. Theoretical mean molecular weight in combustion chamber and at 
nozzle exit for liquid hydrazine with liquid fluorine. Isentropic expansion to alti 
tude indicated assuming equilibrium composition. 
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(d) Combust ion -chamber pressure, 1200 pounds per square inch absolute. 

Figure 6. - Concluded. Theoretical mean molecular weight in combustion chamber and at 
nozzle exit for liquid hydrazine with liquid fluorine- Isentropic expansion to 
altitude indicated assuming equilibrium composition. 


to;- 










Specific heat at constant pressure, c^, cal/(g)(°K) 






O 1.4 [ 


Pres sure , 
. atm 


Stoichiometric 

ratio 


Fuel in propellant, percent "by weight 

(b) Combustion -chamber pressure, 300 pounds per square inch absolute. 

Figure 7. - Continued. Theoretical specific heat at constant pressure of combustion 
products (including energy of dissociation) of liquid hydrazine with liquid fluorine 
, Isentropic expansion to pressures indicated assuming equilibrium composition. 
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(a) Combustion-chamber pressure, 150 pounds per square inch absolute. 


Figure 8. - Theoretical coefficient of viscosity of combustion products 
of liquid hydrazine with liquid fluorine. Isentropic expansion to 
. pressures indicated assuming equilibrium composition. 
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(d) Combustion-chamber pressure, 1200 pounds per square inch absolute. 

Figure 8. - Concluded. Theoretical coefficient of viscosity of com- 
bustion products of liquid hydrazine with liquid fluorine. Isen- 
tropic expansion to pressures indicated assuming equilibrium 
composition. 
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Figure 9. - Theoretical coefficient of the 
tion products of liquid hydrazine with 1 
expansion to pressures indicated assumin 
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(b) Combustion-chamber pressure, 300 pounds per square inch absolute. 


Figure 9. - Continued. Theoretical coefficient of thermal conductivity of combustio 
products of liquid hydrazine vith liquid fluorine. Isentropic expansion to pressu 
indicated assuming equilibrium composition. 
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Figure 11. - Theoretical specific impulse for chamber pressure of 
300 pounds per square inch absolute and exponent n for equation 

1 = *300 (V 300 ) 11 for hydrazine with liquid fluorine. 

Isentropic expansion to expansion ratio indicated assuming equili- 
brium composition. 
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Figure 13. - Theoretical ratio of* nozzle-exit area to throat area for 
chamber pressure of 300 pounds per square inch absolute and exponent 
n for equation S e /S t = ( s e/ s t)300 (P c / 300 ) n for liquid hydrazine 

with liquid fluorine. Isentropic expansion to expansion ratio indi- 
cated assuming equilibrium composition. . 
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.gure 14. - Adiabatic combustion and isen- 
tropic expansion from P c to P e - 
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Figure 16. - Nonisentropic expansion from 
P c to P e . 
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Figure 15. - Loss of heat in combus- 
tion chamber followed by Isentropic 
expansion from P c to P e - 



Entropy, s, cal/(g)(°K) 


Figure 17. - Adiabatic combustion and 
isentropic expansion from P c to 
p x ; removal of heat q. at constant 
p x ; isentropic expansion from P x 
to p e . 
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